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We introduce and discuss a general criterion for the derivative pricing in the general situa- 
tion of incomplete markets, we refer to it as the No Almost Sure Arbitrage Principle. This 
approach is based on the theory of optimal strategy in repeated multiplicative games orig- 
inally introduced by Kelly. As particular cases we obtain the Cox-Ross-Rubinstein and 
Black-Scholes in the complete markets case and the Schweizer and Bouchaud-Sornette 
as a quadratic approximation of our prescription. Technical and numerical aspects for 
the practical option pricing, as large deviation theory approximation and Monte Carlo 
computation are discussed in detail. 

1. Introduction 

The classical ingredients for the derivative pricing and hedging are absence of costs, 
efficient and complete market. According to these hypothesis, the absence of arbi- 
trage opportunities determines the price. 

The moat general formulation of the no-arbitrage argument is due to Harrison 
and KrepsE3. They show that a price system (under certain restrictions such as no 
trading costs) admits no arbitrage opportunities if and only if all price j^^^eiLaie 
martingales. This theory is discussed in several excellent monographs I13'E3ejO'E3. 
Nevertheless the martingale is not unique except for a complete market. In this case 
it is possible, using options and shares of a single stock, to build up a portfolio with 
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deterministic return. Then, the capital growth rate equals the bank interest rata. 
The most successful (and famous) applications are the Black and Scholes formula 13 
for continuous time and the Cox, Ross and Rubinstein for binomial processes c3. 

Unfortunately, real stock markets are not complete and a universally accepted 
pricing procedure is still lacking. Nevertheless, the mean-variance methods seem 
to have the largest consensus BE3. In these theories, the optimal hedge portfolio is 
found by assuming that investors are risk averse. One of the main results is that 
the expected capital invested in this portfolio grows according to the bank interest 
rate. A serious conceptual and practical problem obviously rises: why should a risk 
averse investor, or any rational investor, put the money in such a portfolio? In fact, 
return is not deterministic and it is expected to be merely equal to the fixed bank 
return. 

In the present paper we shall consider an incomplete market under the standard 
assumption of no transaction costs. We state an argument which is able to specify 
the appropriate martingale out of the many possible ones. Considering the deriva- 
tive as a possible asset in the diversification of a portfolio we take the point of view 
of a speculator, who is a market operator interested in the "best" investment in 
the long run limit. We shall show in section 3 in which sense this point of view is 
reasonable. 

The criterion is based on the theory of optimal gambling due to Kelly El and 
hereafter we will refer to it as the Principle of No Almost Sure Arbitrage. Some 
results of this paper have been already presented in (Aurell et al. 98). 

Remark that our methodology differs from the standard approach, where one 
takes the point of view of the investor who tries to minimize risks with an appro- 
priate hedging procedure. In the context of a speculator, we consider in this paper, 
the word "hedging" indicates the portfolio strategy he chooses to reach his aim. 

The paper is organized as follows: in section 2 we discuss the Kelly theory in 
the general framework of incomplete markets. In section 3 we state the Principle 
of No Almost Sure Arbitrage and we obtain the derivative price. In section 4 we 
show the almost sure consequences of an incorrect pricing. Section 5 is devoted to a 
discussion of relations and differences with other pricing procedures. In particular, 
we obtain the Cox- Ross- Rubinstein and Black-Scholes prices in the case of complete 
markets and we derive the Schweizer and Bouchaud-Sornette price as a quadratic 
approximation of our prescription. We discuss the large deviation approach in 
section 6 and we develop a Monte Carlo method in section 7. Finally, in section 8 
we summarize and discuss our results. In the Appendices we discuss in detail some 
technical points: the case of time correlated returns and heavy-tail distributions. 

2. The Kelly theory of Optimal Portfolio 

In this section we summarize Kelly's theory of optimal gambling @. Kelly originally 
was looking for an interpretation of Shannon's Information Theory E3 outside the 
context of communication. Later Breimantm reconsidered Kelly's theory as a model 
for optimal portfolio in a stock market. For a recent review of growth-optimal 
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investment strategies see I 

In the paper we shall consider, except where differently stated, a discrete time 
price movement of stock (or some other security) described by 



u t S t . (2.1) 



where time is discrete, St is the price at time t of a share and the Ut's are indepen- 
dent, identically distributed (i.i.d. hereafter) random variables. 

Let us consider a market operator endowed with a wealth Wo at time zero, who 
decides to gamble on such a stock many times. He invests at each time a fraction 
I of his capital in stock, and the rest in a risk-less security, i.e. a bank account. A 
non-zero risk-less rate corresponds to a discount factor in the definition of the share 
prices, and can be accounted for by a redefinition of the u t 's. We set for simplicity 
the risk-less rate to zero. We will write the result for the general case at the end. 
At time t the market operator will hold a number IWtj St of shares, and his wealth 
at successive instants of time follows a multiplicative random process 

W t +i = (1 - l)W t + lu t W t = (1 + l{u t - 1)) W t . (2.2) 

As a consequence of the large numbers law the exponential growth rate of the wealth 
is in the large time limit a constant with probability one. That is, 

1 

= lim - log -f = E^\og (1 + l(u - 1))] (2.3) 

T^OO T Wq 

holds for almost all realizations of the random variables ut, where 

E?[ (■)]=/ (•) p(u)du . (2.4) 



and p(u) is the probability density for the random factors u. 

The optimal gambling strategy of Kelly consists in maximizing A(/) in (2J3) by 
varying I. We call speculator a market operator who follows this rule. The solution 
is unique because the logarithm is a convex function of its argument. Let us discuss 
which values of I are reasonable in our problem. First, the optimum I must be such 
that l + l(u — 1) is positive on the support of u. Second, we must decide if borrowing 
of cash or short selling of stock is allowed. In the original formulation of Kelly it is 
not, but here it is useful to allow I to take any finite positive or negative value, and 
look for the maximum of X(l). 

The desired strategy is specified by the only finite /* which solves 

and the maximum growth rate is 



A* = E p [\og (1 + l*(u- 1))] 



(2.6) 
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Let us define 

q{u) = l + Z-(u-l) ' (2J) 
where i* is the optimal fraction. It is easy to show that q(u) is a probability density. 



Since l + l*(u — 1) is positive also q(u) is positive. From (2.5) and the normalization 
of p(u) one has 

r q(u)du = J q(u)(l + l*(u- l))du= 1 . (2.8) 

The q(u) is a new probability with respect to the given stock. Furthermore, accord- 
ing to this probability the associated multiplicative process is a martingale. From 
the definition 

E*[ (■)}=[ (■) q(u)du (2.9) 



and using (2.5) it is immediate to obtain 

E q [u] = 1 . (2.10) 



Maximizing any monotonic convex function gives the same formal result as (2.3) 
i.e. a unique solution and, in absence of constraints, the price is a martingale. 
However, Kelly has shown that repetition of the investment many times gives an 
objective meaning to the statement that the growth-optimal strategy is the best, 
regardless to the subjective attitude to risk or other psychological considerations. 

The generalization to a speculator who can diversify his investment on N risk 
assets is straightforward and the interested reader can find it in u. In the next section 
we shall consider the particular case of a two risky assets portfolio composed by a 



share and a derivative written on the share. We shall show that the probability (2.7) 
plays a privileged role for pricing derivatives. 



3. Principle of No Almost Sure Arbitrage 

We consider now the problem of pricing a derivative in incomplete markets. As in 
the previous section the speculator can invest his capital only at discrete times. 

The speculator's portfolio is composed at time t — by a risk-less security, a 
stock and a derivative written on the same stock with maturity time T and strike 
price K = kSo. 

At time T the derivative expires and the speculator decides to buy a new deriva- 
tive with the same k and maturity IT . Let us focus on the investment at time t in 
order to fix the notation. We call C t the value of derivative at time t and C t +i(u t ) 
the value at the following instant of time t + 1. Finally 

ft(ut) = ^ (3.11) 

is the return of the derivative. The notation ft(ut) is to stress that we are dealing 
with a derivative, i.e. the return of this asset depends on the return of the share u t . 



A general methodology to price and hedge derivatives in incomplete markets 5 



The simplest example is an European call option with a strike price kSo, one time 
step before the expiration date i.e. t = T — 1; in this case one has Ct(ut — 1) = 
{ut — 1St-i — kSo} + . 

At time t the speculator diversifies his portfolio by investing a fraction l t of his 
wealth in shares and a fraction dt in derivatives. If the speculator keeps some money 
Ws in stock and Wd in derivatives, after the market fall out the two capitals are 
worth respectively WsUt and Wnftiut) where the same random variable u t appears 
in both expressions. The total wealth is therefore changed to 

W t+1 = (1 + l t (u t - 1) + d t (ft(u t ) - 1)) W t ■ (3.12) 

We shall consider hereafter, as in section 2, the return u± an i.i.d. random 
variable. This assumption is relaxed in appendix A where the markovian case is 
treated. 

The return of the derivative ft is a function ofT — tifO<<<T and of the 
random variable Ut- Let us stress that if Ut is an i.i.d. random variable the returns 
ft are periodic of period T. Hence the fractions {h}t and {dt]t take on only T 
different values. 

We consider the case of a speculator who plays the game every time and repeats 
his investment many times. In this situation his wealth will almost surely grow at 
an exponential rate which is 

Kihh, {dt}tl {fth) = lim - V log (1 + l n (u n - 1) + 4(/»K) - 1)) = 

r — >oo T — 
n=0 

T-l 

= - E p [log (I + l t (u t - I) + dt(ft(u t )-l))] .(3.13) 



t=o 



The speculator is then interested in the fractions attaining the optimal growth rate 
of the capital ({ij }*=o...t-i, {<^t}t=o...T-i)- This set of values is unique because of 
the convexity of the logarithm given the set of the derivative returns {/t}i=o...r-i- 

We are now ready to state the Principle of No Almost Sure Arbitrage (hereafter 
NASA) introduced in la: prices of derivative securities must be such that a spec- 
ulator cannot construct a portfolio out of combinations of the derivatives and the 
underlying security which grows almost surely at a faster exponential rate than a 
portfolio containing only the underlying security. 

In other words, the NASA Principle says that the returns {ft}t must be such 
that the global maximum of ( |3.13] ) equals A*, i.e. 

A({i?}*,to}t;{/t},)=A* (3.14) 



where A* is the Kelly optimal rate of equation (2.6). Since one trivially has \({lt = 
l*}t, {d t = 0} t ; {ft}t) = A*, the maximum must be in (1% = l*,d^ = 0), because of 
the uniqueness of the maximum of a convex function. 

Let us give a simple intuitive interpretation of this fact. If dl is larger than 
zero all speculators would like to buy the derivative in order to achieve a larger 
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exponential rate of their capital. As a consequence the derivative price rises and 
its return rate falls. The resulting fraction djf decreases. On the other hand, if 
is less than zero, then all speculators want to go short of derivatives causing their 
price to fall and the return rate to rise so that djT tends to grow. 

From a technical point of view, the idea of setting d£ to zero, is closely re- 
latedto the Samuelson's methodology E3 of warrant pricing in the 'incipient case', 
see tEl[chapter 7]. We apply such a principle here to derivatives and not only to 
warrants. Furthermore, the condition d* t = stems from a general principle. The 
incipient case is intrinsically important and not only a convenient assumption to 
simplify the computations. 

The condition for is 

= Wt \ {lt ^ }u{dt =o }t = E p l 1 + Hu _ 1) ] ■ (3-15) 



This equation is identical to ( |2.5| ) of the Kelly model and we find again the proba- 
bility q ( |2.7|) . The equation specifying dj? is 

n dX({lt}t,{dt}f,{ft}t) l FPr (ft(u t )-l) 

° = dd t k^h,{^o }t = m i + nut _ 1} } ■ (3-16) 



Then using (^7J) 

C t =E q [C t+ i(u)} , (3.17) 

which states that the NASA price for derivatives is a martingale with respect to q: 
in the following we refer to it as the equilibrium price. In appendix B we deal with 
some simple cases where the q can be computed analytically, and we show that even 
a p(u) very wildly divergent for large u will lead to a q(u), which has at least finite 
first and second moments. 

We have derived the rule which gives the price at time t given the price a time 
t + 1. In particular, the option price at time T — 1 is 

Ct-i(St-i) = E«[C T (uS T -i)] ■ (3.18) 

since the price of the derivative is known at the expiration time T. If we iterate 
backwards the above equation we obtain 

Ct(S t )=E Q [C T (U T - t S t )} (3.19) 

where the expectation value [■] is taken with respect to the random variable 

a T 

u T - t =-f= n ^ ( 3 - 2 °) 

^ s=t+l 



which is the product of T — t independent identically g-distributed variables u s . 
The probability distribution Qt-^StIS*) of the product can be obtained as the 
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convolution of T — t probability distributions of a single variable. The price depends 
on the actual share price St at time t. In the particular case 

C t {St) = {St - kS } + (3.21) 

we obtain 

C t (S t )=EQ [{U T ^ t S t - kS a } + ] , (3.22) 

which is the NASA price of an European call option. 

We have taken the risk-less rate to be zero. A non-zero risk-less rate can be 
reintroduced as a discount factor in the share prices. Actually, having defined u 
as u/r we can reproduce all the computations for the new stochastic variable. We 
start from the effective probability of the u and we find the probability q(u). Then, 
for r > 1 the European call option price is 



C t (S t ) =E<3 



Ur-tSt IfT—r 



(3.23) 



where U T -t = IlLt+i 

Before we turn to the consequences of a derivative price different from the one 
imposed by the NASA Principle, let us pause and comment. We are interested 
in the derivative value assuming that the speculator knows but cannot change the 
distribution of stock price returns. This is of course an idealization but it looks a 
reasonable starting point. We stress once again that the Principle, as it is stated 
above, can be applied only to an asset and a derivative written on it, because both 
of them use the same source of information: the distribution of the returns of the 
underlying. Surely it cannot be applied to determine the stock price. In this case 
two different financial objects are involved, a risk-free interest rate and a risky 
asset. It is not reasonable to have the same growth rate in both cases, otherwise 
nobody would invest in the second one. However we expect to observe in liquid 
markets, the tendency of the average discounted return to become closer and closer 
to zero. This is true for example in the case of currency markets. In fact after 
1989 for the currencies of the most developed countries the purchasing power parity 
is well verifiedji.e. the exchange rate of two countries should adjust according to 
relative prices □ A null average return, for periods of one or few years, is a simple 
consequence of the purchasing power parity and it implies that the optimal fraction 
in a stationary policy is zero. 



4. Arbitrage in non equilibrium case 

We show here that an incorrect pricing of a derivative allows for almost sure ar- 
bitrage i.e. a rate of capital grow larger than the one obtained with the optimal 
portfolio of shares and risk-less securities. 

*The observed deviation from this behavior is generally due to macroeconomic reasons, such as 
the rapid depreciation of Italian Lira in the wake of ERM crisis in 1992 or the large appreciation 
of UK exchange rate in the eighties with the development of the North Sea oil fields £3. 
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Let us suppose that the derivative is not correctly priced (i.e. it is not at its 
equilibrium value C t *): 



C t = C* + AC t with t = 0, T - 1 



(4.24) 



Given the price of the derivative the problem is to find the optimal growth rate 
of the capital for a mixed portfolio composed by a risk-less security, a stock and 
the derivative written on it. Assuming that AC t is small and expanding the growth 
rate up to the second order in Sit = h — I*, 8d t = d t and ACt we find that the 
optimal portfolio corresponds to 81% and 8d% given by 



where the matrix is 

EP 
E p 



si: 



Sd* t = 



acv r 



(*) 



Ct detTW 



Ac t r 



ii 



c* detrW 



r (t) = 

The corresponding maximal rate is 



(1 + J.( U _1))2 

(ft(u t )-l)(u t -l) 



EP 
E p 



(/ t («t)-i)(«t-i) 



A* 



1 T_1 



(i+;«(« t -i))' 

(/.(«.)-i) 2 
(i+;«(« t -i))^ 



,(t) 



2T f-' V Ct / detrW 



(4.25) 
(4.26) 

(4.27) 
(4.28) 



Remark that and detr'*^ are positive quantities. Therefore, the correction to 
the equilibrium optimal growth rate A* is always positive both for under- valued and 
over- valued derivatives. 



5. Comparisons with other approaches and limiting cases 

In this section we compare our approach with previous pricing procedures both in 
case of complete and incomplete markets. 

5.1. Classical pricing prescriptions for complete markets 

The above proposed pricing procedure agrees with no-arbitrage pricing in the case 
of complete markets. It is, therefore, instructive to carry through the calculations 
for the dichotomic Cox-Ross-Rubinstein case and the Black-Scholes continuous time 
limit. 

Following Cox-Ross-Rubinstein we assume that the share price can go up by a 
factor u u with probability p and down by a factor Ud with probability (1 —p). We 
can safely assume that the risk-free interest rate r equals unity. The case r > 1 can 
always be recovered by simply replacing u with u — u/r. The probability density 
is therefore: 

p(u) =pS(u- u u ) + (1 - p)8{u - u d ) (5.29) 
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where the £(•) are Dirac delta functions. 

The optimization equation is then solved as 



r = _p_ _ (1-pl 
1—Ud u u — l 



which leads to a probability density q(u) having exactly the same form as ( 5.29| ). 
the only difference being that the probability p is replaced by q: 

v 1 — U r l ._ 

(5.31) 



l + l*(u u -l) U u —Ud 

which is the Cox-Ross-Rubinstein probability q. The new probability q is indepen- 
dent of p: a result which suggests that this case is somehow atypical. 

Consider now the continuous limit case of Black and Scholes. The prices of shares 
are continuously monitored and the financial operator is allowed for a continuous 
time hedging. In an infinitesimal time interval dt the share changes by a factor 

u = exp{-qdt + Adw} (5.32) 

where rj and A are constants while dw is a random increment with vanishing average 
and variance E p [(dw) 2 ] = dt. This infinitesimal time expansion is equivalent to say 
that 

E p [logu] = rjdt , Var p [log u] = A 2 dt (5.33) 

while it is assumed that higher moments of log(it) have expected value of higher 
order in dt. 

With the above assumptions one can easily solve the equation for I*, in fact, up 
to order dt equation for I* is 

E p [{u I* (u-1))} =0 (5.34) 

which gives 



/ =2 + 1- ( .-,. 8 r„ 



Using ( [5.35 ) it is easy to verify that 

A 2 

E q [logu} = — —dt ; Var q [\ogu] = A 2 dt (5.36) 

and therefore for a finite time T the variable Ut — St /So is distributed according 
to a lognormal martingale which follows the same prescription of Black and Scholes. 

Notice that the continuous case of Black and Scholes refers, as well the Cox- 
Ross-Rubinstein case, to a complete market. In both cases it is possible a perfect 
hedging, i.e. a complete replication of the option price process by a combination 
of shares and risk-less investments. This fact is not surprising since the continuous 
gaussian case can always be recovered from the dichotomic case, when the the time 
step becomes infinitesimal. 
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It is useful to stress that, even if the price with our principle is the same of the 
classical approaches in complete markets, this does not mean that the strategy is 
the same. The speculator, in our case, does not worry to hedge continuously in 
order to obtain a risk-less portfolio but he follows the strategy of trying to do the 
best with the information he has i.e. he attempts to get the largest almost sure rate 
of exponential growth. Even in the complete market case, the speculator obtains 
an almost sure exponential growth rate A generally larger than the risk-less interest 
rate r (obtained by the classical pricing prescriptions), explaining why one should 
invest in a portfolio which needs an active and frequent trade. 



5.2. Schweizer and Bouchaud-Sornette approach as quadratic approxi- 
mation 

In the Schweizer and Bouchaud-Sornette i approach the option price is ob- 
tained searching an hedging strategy which minimizes a quadratic risk function. 
This is a natural generalization of the Black and Scholes idea of looking for a zero 
risk strategy (perfect hedging) to the case in which this is not anymore possible (in- 
complete markets). In the limit of continuous time the usual Black-Scholes theory 
is recovered. The problem is that the above approach gives negative values of the 
optkm price for large St I23eH as shown in a different context by Dybvig and Inger- 
soll 113. Nevertheless, negative prices appear in very unrealistic situations while, for 
more realistic cases, the formula gives sensible values. 

We show in this section that the Schweizer and Bouchaud-Sornette price can be 
considered a quadratic approximation of the NASA one. 

Let us first recall their result in a convenient form. It has been shown in lid 
that their price can be expressed in the same way of equation (3.19) where the 
expectation is now constructed using the pseudo-martingale probability 

q(u)=p(u){l--^[(u-l)-u]} (5.37) 

where u = E p [u — 1], a 2 = Var p [u — 1]. 

According to a rather general theory, one can show llil that the martingale prob- 
ability associated to an utility function U(W) in a two times problem is 

g ( u ) = p (u) wE^tll . (5.38) 
E p [^U{W)U.] 

where 

W= (l + i(u-l))W Q (5.39) 

is the capital after the investment and Wq before it; I* is the optimal fraction 
associated to the utility U. 

From equation ( 5.38 ) is clear why negative prices appear: this happens only 
when the utility is a decreasing function of the capital. 

Let us introduce a quadratic utility function of the capital W 

U quad (W) = aW ~~{W -Wf (5.40) 



A general methodology to price and hedge derivatives in incomplete markets 11 



where W is the expected capital and a and (3 two positive parameters. 
Inserting this utility JJ quad in equation ( J5.38D one gets 

q{u)=p{u){l-l*W -[{u-l)-n}} . (5.41) 
a 

The optimal I* can be derived simply by the imposition of the martingale property 
or equivalently from the equation 

rw - = 4 . ( 5 - 42 ) 



obtaining that the above pseudo-probability coincides with ( 5.37 ). 

Notice that U quad (W) is nothing but a Taylor expansion up to the second order 
of a generic utility function around W and then also of the logarithmic case of the 
speculator we have considered. Our price and the Schweizer and Bouchaud-Sornette 
one coincide for all practical purposes in mild periods but can lead to different results 
if strong fluctuations of the returns are present (i.e. near financial crisis). 

6. Large deviations 

In section 3 we have shown that the NASA price of a derivative is determined by a 
probability distribution Qt{St\Sq) which is related to the one-time step equivalent 



martingale probability q defined in (2.7) 



The probability Qt(St\Sq) is constructed by compounding the one time step 
probability q. If u can take only a finite number of values N, then Qt{St\Sq) 
can be written as a multinomial formula, a simple generalization of the Cox-Ross- 
Rubinstein case with binomial coefficients. 

Unfortunately, if the number of time steps is large, it is not possible to compute 
explicitly the multinomial average which gives the price because it involves 0(N T ) 
operations, and some other approach must be found. Therefore we need some 
analytical and numerical techniques for an estimation of the price. In this section 
we discuss the large deviation approach, next section is devoted to the Monte Carlo 
method. 

The natural answer to the problem of a correct approximation of the price comes 
from large deviation theory. To state the problem we have to define the proper large 
deviation variable 

z = ^logS T /S Q (6.43) 

whose probability distribution Qt(z) is trivially obtained from Qt(St\Sq) as 

Q T {z) = T exp{zT} Q T {exp{zT}) . (6.44) 

In a nutshell, the mathematical essence of the large deviations theory @ is the 
existence ( under suitable hypothesis: essentially the finiteness of all moments of 
exp(Tz)) of a convex function G{z) (usually called Cramers function) such that 

Qt(z) = $ r (a0 exp{-G(z)T} (6.45) 
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where $t(z) is subexponential for large T. 

It is often reported in literature that the distributions are heavy-tailed, and 
therefore it could be doubtful the existence of Cramer function because the hypoth- 
esis is not verified. However distributions obtained by historical sequences always 
have all finite moments, this is a consequence of the presence of a natural cut-off 
which can also be reflected at the levaLpf models c3. 

Because of the Oseledec theorem one has that z, in the limit of very large 
T, is almost surely equal to its expected value A, therefore G(z) has to be strictly 
positive except for z = A where it vanishes. 

In order to find out in practice the Cramer function it is sufficient to compute 
the scaling exponents L(n) for the moments of St 

L(n) = lim hog [(S T /S ) n }. (6.46) 

T — >oc 1 

where 

E Q [(S T /S ) n } = J Q T (z)exp{zTn}dz . (6.47) 

in fact, the Cramers function and the L(n) are related via a Legendre transforma- 
tion: 



G(z) = max n [nz — L(n)] = zn*(z) — L(n*(z)) (6.48) 

where n*(z) is the value where the maximum is realized. 

The lognormal approximation can be recovered keeping only the first two terms 
of the Taylor expansion of L(n) around zero. One obtains 

G(z) = ^(z-\f (6.49) 

where 

A ee E q [logu] , A 2 = Var q [log u] (6.50) 

which is a good approximation only for relatively small fluctuation of St around 
its typical value exp AT, In spite of its relevance in many fields the lognormal ap- 
proximation is from a mathematical point of view rather peculiar. This is basically 
due to the fact its moments grow too fast and therefore the Carleman criterion does 
not hold A simple way is to consider the correction to the parabolic shape of the 
lognormal approximation is to consider the Taylor expansion of order R of L(n) 

R 

L{n)^Y. X o nl ( 6 - 51 ) 

3=1 

that implies 

R. 

J=2 
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where 

.93 = -h 

nr - 1 ( A 3 A ^ A 3 J- X A 
■95 - Af ^AJ — + Xj 

_ 1 / 7 A* 7 AlA 4 1 /A* , , \ A 6 \ 
■96 - Af ^3 Af " 2 "AT + Aj ^X + A 5 A 3 J - — J 



The above expansion of £(n) around n = is quite reasonable because in this 
way one has a rather good approximation of the probability distribution around the 
maximum. On the other hand, at least in non trivial cases, one does not obtain good 
results for the option pricing, even if probability normalization factor are properly 
taken into account (see figure ||). 

In order to understand this fact it is enlightening to compute the scaling expo- 
nents L(n) using ( |6.45|) and ( |6.46 ). 

EQ[(S T /S ) n ] ~ exp{T(nz» - G(z*(n)))} (6.54) 

where the symbol ~ indicates logarithmic equivalence and z*(n) is the maximum of 
nz — G{z). The most important contribution to this integral comes from rare events 
of exponentially small probability such that z — z*(n) ^ A. Only when n — the 
integral is dominated by the most probable events z*(0) = X. 

Now we can easily see how the above considerations enters into the problem of 
pricing options. 

Let us notice that the price Co of an European call option can be written as 

Co = So J dzQ T {z) {expTz - k} + = S h - KI 2 (6.55) 

where 

h = Ix\ osk dzQ T {z) expTz I 2 = J± logk dz QT(z) . (6.56) 

Then the origin of the problems both for the naive Monte Carlo (see next section) 
and large deviation becomes clear: the most important contribution to the integral 
I\ is given by the exponentially rare events for which z = z*(l). It is then important 
to observe that because of the martingale property 

Clr(z) = QT(z)expTz (6.57) 

is a probability distribution. Since the most important contributions to I\ come 
from the most probable events of the Qt(z) distribution, we can simply repeat the 
large deviations expansion of L(n) around n = for this new probability distribution 
while li is computed in the previous way. 

It easy to check that L(n) computed with respect to iV exactly equals L(n+1) 
computed with respect to Qt- Therefore the expansion with respect to n = for 
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Fig. 1. Option prices C rescaled by the initial value of the stock So as function of log(K / So) 
for T = 30. The returns of the price over one elementary time step take three discrete values 0.2, 
1 and 2.5 with probability 0.15, 0.15 and 0.7 respectively. We compare the exact prices, and the 
Black-Scholes ones with the large deviation prescription with an expansion of L(n) around with 
R = 2 (lognormal approximation) and R = 6 (see text). 
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Fig. 2. Option prices C rescaled by the initial value of the stock So vs. log(K/So) for T = 30. 
The returns of the price over one elementary time step are the same of figure 1. We compare the 
exact prices, and the Black-Scholes ones with the large deviation prescription described in the text 
with R = 2, 6. 
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Fig. 3. Option prices C rescaled by the initial value of the stock So as function of \og(K / Sq) 
where K is the strike price, T is chosen equal to 30. The returns of the price over one elementary 
time step are the same of figure 1. We compare the exact prices, the ones obtained with the 
naive Monte Carlo algorithm with 1 million trials and the Black-Scholes prescription. The values 
have been intentionally chosen unrealistic in order to stress the differences between the different 
approaches. 

the new probability distribution corresponds to an expansion around n — 1 for the 
old one. 

All this integrals are not very sensible to the choice of the sub-exponential func- 
tion <&t{z) which can be fixed by means of a normalization constant or, if more 
computational precision is needed, as 



where is the 2 nd derivative of the function G. 

In figure || we show the results obtained in this way for the cases with R = 2 
(lognormal approximation) and R = 6. We observe a quite good agreement with 
the exact prices. 

As a final remark we want to stress that the lognormal approximation does 
not coincide with the Black and Scholes except in the continuous time limit. The 
differences can be appreciated in figure [l] and ||. 

7. Monte Carlo 




(6.58) 
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The idea to use the Monte Carlo method c3, to compute quantities which de- 
pend on markovian processes, is very common in both the physicist and economic 
community Instead of computing an average of a function A(St) according the 
probability Qt(St\Sq): 

EQ[A(S t )} = I A(S T )Q T (S T \S Q )dS T , (7.59) 

one performs a large number M of trials (i.e. M different realizations of the random 
process St according to the probability Qt) and then compute the average as 



4 M) = ^E^(4 j) ) (7-60) 



3=1 

where S^ ' is the j — th realization and E®[A(St)] = limjw-^oo Ajf 1 '. We discuss 
the case in which the u can take only a finite number of values u^', i = 1, 2, N 
with probability pi . The j — th realization is obtained in the following way. From 
So one constructs S [ f) as S[ f) = S ui, where u\ is drawn with probability ( |2.7| ), 
one extracts in the same way u-i and so on; with this procedure one obtains S\f* , 
■ ■ ■ and finally Sip . In this way the number of operations involved is O(MT), 
while the error estimation using standard variance considerations is 



a4 m > 



\ 



M 

_L-2(A(^)-4^)a (7.61) 

j'=i 

For example in figure ^ we show the price of an European call option as a 
function of a trinomial probability distribution. The prices obtained with Monte 
Carlo simulations with 1 million trials are compared with the Black-Scholes ones 



and the exact ones given by Eq. (3.22) 



The large errors involved in this approximations are mainly due to the fact that 



the most important contribution in (3.22) comes from the tails of the distribution 



Qt{St\Sq). We need then a good control of the tails; this is provided by the large 
deviations theory or a suitable Monte Carlo. 

Taking into account the results of the previous section we have now all the 
ingredients to build a reasonable Monte Carlo algorithm. As previously shown 



in (6.55) the price is a linear combination of two integrals, I\ and li- The most 
important contribution to I± come from the most probable events of fix- Let us 
write 

dzClr(z)= n T(S T /S a ) (7.62) 

r lo g fc S T >K 

where fi, in the case u assumes only a finite number of values N, is the convolution 
of a number T the one-time probability defined by 



w(«(*') = u (i )<z(u«) 



(7.63) 
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Fig. 4. Option prices C rescaled by the initial value of the stock So vs. log(K / Sq) for T = 30. 
The returns can take the values uv n ) = u^a n , with equal to 0.4, a equal to 1.2 and the 
index n ranging from zero to TV = 30. The probability p(«( n )) is Cjv(u (n) )-' 3 , with fi = 0.5 
and Cjv a normalization constant. We compare the large deviations (R = 6) and Monte Carlo 
approximations and the Black-Scholes prices. The Monte Carlo simulation has been performed 
with 10 7 trials. 
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Fig. 5. We compare the same quantities of figure |4| with the exact prices in the case of the 
trinomial returns of figure 1. The Monte Carlo simulation has been performed with 10 5 trials. 
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that is properly normalized because of the martingale property. The algorithm 
follows then the same steps of the simple Monte Carlo, using to compute the integral 
I\ the w-probability instead of the o-probability. To improve the results we have also 
used the control variate technique u. Even in the extreme situation of a very large 
Var[u] such as in the case of a truncated Levy distribution of figure ^ we obtain 
good results both for the Monte Carlo and the large deviations approximation. In 
figure H we show the results obtained the trinomial distribution of figure 1. 

8. Final Remarks 

In this paper we have discussed a general criterion to price derivatives. This prin- 
ciple (No Almost Sure Arbitrage) basically says that it is not possible to build a 
portfolio including derivatives which grows almost surely with an exponential rate 
larger than one with only the underlying securities. 

This approach holds for the general case of incomplete markets and for a generic 
probability distribution of the return. Let us stress that it is possible to apply the 
procedure also for returns correlated in time, e.g. a Markov process. 

In the cases of dichotomic distribution for the return and continuous time limit, 
our method gives the same results of Cox-Ross-Rubinstein and Black-Scholes re- 
spectively. If p(u) is not too spread, the results obtained with our approach are 
rather close to those ones of the Schweizer and Bouchaud-Sornette theory. This 
is not true in presence of strong fluctuations for the return. In this case an accu- 
rate computation of the option price can be obtained in the framework of the large 
deviations theory which systematically takes into account the deviation from the 
lognormality or via a properly realized Monte Carlo computation. 
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Appendix A 

In this appendix we consider a Markovian model for the price movement of the 
stock (i.e. equation ([O])) in order to extend our price prescription to this case. 
Now the u t are random variables with a Markovian probabilistic rule. For sake of 
simplicity we discuss the case of Markov chain, i.e. u t can take only N different 
values it' 1 ), u^ 2 \ ■ ■ ■ u^ N \ 

Let us remind that a Markov chain is completely characterized by its transition 
matrix 
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Pi^j = Prob(u t+1 = u {]) \u t = u {l) ) . (A.l) 

In a realistic market it is natural to assume that the Markov chain is ergodic and 
therefore there exist an invariant probability: 

N 

Pj = ^PiPi^j ■ (A.2) 
For ergodic chain the convergence to the invariant probability is exponentially fast 

(P n )^3 = Prob(u n = u {l) K = = Pj + 0(e- an ) (A.3) 

where a is given by the second eigenvalue of the matrix Pi—>j ■ 

The optimization strategy for the shares problem is rather obvious. If at time t — 
1, Ut-i = then the probability to have u± = uW is nothing but Pi->j. Therefore 
the optimal fraction of capital in stock is obtained from the maximization of 

JV 

Y,Pi^±og[l + k(uW-l)}. (A.4) 
i=i 

Of course U can change at any i. In Kelly's paper is discussed the relation between 
the Shannon entropy of the Markov chain and the optimal exponential growth rate 
of the capital. 

Following the same step of Kelly's strategy described in section 2 it is straight- 
forward to construct from Pj_>j an equivalent martingale with transition probability 

Qi — >j • 

P — (A.5) 



1 + ' 



where I* is given by the maximization of the quantity (A.4). It is easy to realize 
that 

E Q {u\u^) = 1 . (A.6) 

If the Markov chain with transition matrix Pi—>j is ergodic then also the Markov 
chain with transaction matrix Qi^j has to be ergodic and one can introduce an 
invariant probability 

N 
i=l 

Also for the new Markov chain the convergence to the invariant probability is ex- 
ponentially fast 

(Q n )i^j = qj + 0(e'^) , (A.8) 
nevertheless, in general, (3 ^ a. 
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Now we have all the ingredients for the construction of QT(St\Sa) since we can 
assign the probability to each trajectory. Then, by repetition of the same steps of 
section 3, one can construct the price of the derivatives. The only difference is that 
now the expectation is taken with respect to the markovian martingale. 

Let us notice that for the dichotomic case N = 2, the present construction gives 
the same probabilistic rule of Cox-Ross-Rubinstein. In this case the Qi^j do not 
depend on Pj_>j but only on the values vP~' and u^ 2 \ The Monte Carlo computation 
is rather simple. In each trial, at time t, if Ut—i — vft', one chooses with probability 
Qi^j the j-th branch of the tree. 

Also the lognormal approximation for the Markovian case is straightforward but 



interesting. Equation (3.4S) holds, now one has 



and 

T 



\= lim ^^[Vlogut] = E Q [\ogu] (A.9) 

t=\ 

A 2 = lim ii^[(f>g Ut - A)) 2 ] . (A.IO) 



T^oo T 

After having defined 

R n ee ^(logu^ - A) (log - X)(Q n )i^j , (A.ll) 

id 

it is a matter of computation to check, using the exponential convergence of R n due 
to the Markovian nature of the process, that 

oo 

A 2 = E Q [(log u - A) 2 ] + 2 ^ i?„ . (A.12) 

n=l 

Let us note the fact that, at variance with A which depends only on the invariant 



probability, for A 2 are also relevant the time correlations. Equation (A.12) shows 
that the "effective" volatility takes into account the time correlation. 



Appendix B 

It is useful to have some cases in which the probability q is easy to handle analyt- 
ically. In this appendix we consider three particular distributions of the return of 
the asset, in the case of small excess rate of return fi = E p [u — 1]. 

Compact support distribution 
We first consider a probability distribution p(u) in a compact support [0, U max \- 
Equation fl2.5p can be rewritten as 

oo 

E p [(-l) n (l*) n {u - 1)™ +1 ] = . (B.l) 

n=0 

We can solve for I* as a power series in the excess rate of return fx 
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r = 4^ + 4m 2 + °(m 3 ) . (b.2) 

<7 Z <J° 

where a 1 and x are the second and the third cumulant of the rescaled return u — 1 
respectively. From equations (2.7) and (B.2) we get the probability 

q(u) = p(u) (l - JL(u - 1) + /i 2 - 1) + ^(« - l) 2 ) + O (m 3 )) , (B.3) 

hence up to order /i the result of the NASA principle is identical to the minimal 
risk prescription of Schweizer and Bouchaud-Sornette. 

Expansion (B.l) is convergent if the optimal value I* is small enough. The 
individual terms in the expansion however exist under the much weaker conditions 
that the respective moments are finite, but the support of p(u) can be unbounded. 

Levy distribution 

There are some known experimental ev idenc es that the one-time return u is well 
described by a proper Levy distribution E3o, but it is still an open problem if the 
variance is finite or not. In this case expansion (B.l) is meaningless. We assume 
here that the returns are described by a probability distribution with a power law 
decay at large arguments 

p(u) ~ fivr 1 "" for u » 1 . (B.4) 

The exponent of the power law v is in the interval (1,2] and j3 is a constant. Hence 
the expected value of u exist and is finite, but all the higher moments are infinite. 

Let us notice that the maximizing I* cannot be less than zero, because the 
support of p is unbounded. Hence, I = is at the boundary of the domain for which 
X(l) is well-defined function of I. 

Let us now assume that I* is small, but not necessarily of the same order as 
in this case the probability distribution q(u) is well approximated by 

q{ ) X ikpM a « » £ ' (R5) 



and then its second moment exists. Using the approximation ( p3.5| ) for the proba- 
bility q and the definition of fj,, after some algebra one obtains 

r /mK^i)\^ (B6) 





In the limit when v tends to 2, I* scales linearly with fj,, in agreement with ( |B.2] ). 
When v is in the interval between 1 and 2, the optimal fraction however goes to 
zero faster than 

Log-Levy distribution 
We construct here the probability q from a p of the form 



^logu)- 1 -" for u»l 
Cu -1 1 logul -1 -^ -5 for u«l 



Pto~< X _i 7 i-i-vC-) <" , (B-7) 
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We will refer to these laws as "log-Levy" laws, since they are distributed as Levy laws 
with exponents and i/ — ) in the logarithmic variable. In this case the expected 
value of the logarithm of the return is finite but the variance is infinite. As the Levy 
laws are stable under addition, so the log-Levy laws are stable under multiplication, 
if the exponents of the power laws, and are in the interval [1,2]. The 

limit case of an exponent equal to 2 corresponds to the lognormal distribution. 

However, this does not mean that the log-Levy laws are necessarily as natural as 
the Levy laws. The divergence at large u in (B.7) implies that the expectation value 
with respect to p of the random variable u does not exist. Therefore the excess rate 
of return fj, is undefined. 

If the probability distribution p is nicely behaved for small u we have that the 
growth rate of the wealth (2.3) exists and is finite in a region around I = 0. We can 
therefore find the maximizing I* , attained at some positive I*. Using the variational 



equation, we obtain the probability (2.7) that can be approximated by ( |B.5| ). Just 
as in the Levy distribution, we have that E q [u 2 ] is finite, and it follows the much 
weaker condition that E q [(\ogu) 2 ] is convergent at large u. 
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